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Abstract 

CePcDAD crystallizes in the CaBe2Ge2-type tetragonal structure ( PA/nmm , 129 ) and under- 
goes a phase transition to the orthorhombic Cmme structure at around 13 K. Its inelastic neutron 
spectra reveal an additional magnetic excitation that was ascribed to electron-phonon interaction 
leading to a formation of a new quantum quasi-bound vibron statě. We present the first-principles 
calculations of the crystal field excitations and lattice dynamics calculations of the phonon dis¬ 
persions to compare with the experimental data. The calculated crystal ňeld energy splitting in 
CePd2Al2 agree well with the model ušed to describe the experimental neutron scattering spectra. 
The hrst excited crystal held level moves to higher energies when undergoing the transformation 
frorn tetragonal to orthorhombic structure, in agreement with experiment. The lattice dynamics 
calculations show that in both tetragonal and orthorhombic structures there are no imaginary 
modes for any q-wave vector within the Brillouin zone and therefore the lattice structures are 
stable. The phonon dispersions and density of states are calculated for both crystal structures of 
CePd2Al2 and its nonmagnetic counterpart LaPd2Al2. The results generally agree well with the 
experimental data including the high phonon density of states around 12 meV. The phonon density 
of states is also ušed to calculate the lattice heat capacity and compared with the experiment. 

PACS numbers: 


1 



I. INTRODUCTION 


The tetragonal CeT 2 X 2 compounds (where R = rare earth, T = transition metal d- 
clcment and X = p-metal) form a large group of materials which exhibit various inter- 
esting physical properties as heavy-fermion behavior, valence-fluctuations, non-Fermi-liquid 
behavior or unconventional superconductivity. 1-5 These interesting properties stem frorn 
competition between RKKY and Kondo interactions and a strong influence of crystal field 
snrrounding the Ce ions. The observation of an additional magnetic excitation in the neu¬ 
tron scattering spectra of CePd 2 Al 2 is another interesting phenomenon. 6,7 The existence of 
three excitations frorn the ground statě of Ce 3+ ions cannot be explained in terms of pure 
crystal field splitting and was ascribed to electron-phonon interaction leading to a formation 
of a new quantum quasi-bound vibron statě similar to CeAl 2 8 ’ 9 and CeCuAR 10 . This strong 
magneto-elastic coupling between orbital and lattice degrees of freedom and the resulting 
magneto-phonon mode háve been described within the Thalmeier-Fulde model in the čase 
of cubic CeAl 2 . n 

CePd 2 Al 2 crystallizes in the CaBe 2 Ge 2 -type tetragonal structure ( PA/nmm , 129) and 
undergoes a phase transition to the orthorhombic Cmme structure at around T s ^. d2Al2 = 
13 K 6,12-14 . Similar structural phase transition is observed also in LaPd 2 Al 2 6,12 excluding a 
magneto-elastic origin of the phase transition. The structural instability is not unusual in 
this family of compounds. Phase transition frorn tetragonal to a lower symmetry structures 
are observed e.g. in CeNi 2 Sn 2 , 15 CePt 2 Sn 2 , 15 CeRh 2 Sb 2 , 16 CePt 2 Ge 2 17 or CePd 2 Ga 2 . 18,19 The 
magnetic properties of CePd 2 Al 2 are characterized by an antiferromagnetic order below Tn = 
2.7 K 13 . The ground statě magnetic structure is an amplitudě modulated wave described by 
an incommensurate propagation vector k = (ó x , \ + S y , 0) with S x = 0.06 and S y = 0.04. The 
magnetic moments order antiferromagnetically within the planeš stacked along the c-axis 
and are arranged along the direction close to the orthorhombic a-axis 7 . The recent neutron 
scattering experiment provided a clear link between the observed magnetic excitations and 
structural properties of CePd 2 Al 2 . The comparison of magnetic excitations measured above 
and below 13 K, where the structural transition frorn tetragonal to orthorhombic structure 
occurs, shows a strong shift of the energy of the hrst excited level to higher energies when 
cooling below 13 K.' The other two excitations are almost unaffected by this transition. 
These results motivated our ab-initio calculations of the phonon spectra and crystal field 
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excitations in both structural phases of CePd 2 Al 2 to explore origin of the experimentally 
observed changes in energy spectra. 


II. COMPUTATIONAL DETAILS 

We performed first-principles electronic structure calculations using the generál potential 
linearised augmented plane wave rnethod (APW + lo, WIEN2k) 20 . The Kohn-Sham equa- 
tions were solved within the generalized gradient approximation (GGA) 21 . The relativistic 
effects were treated using scalar relativistic approximation. Atornic sphere radii (AS) of 
2.5, 2.1 and 1.8 Bohr rádius (1 Bohr rádius = 52.9177 pm) were chosen for Ce, Pd and AI, 
respectively. About 1400 linearized augmented plane waves (140 per atom) were ušed in the 
interstitial region and the highest value of 12 in the expansion of radia! wave functions inside 
the AS to represent the valence states. The correlated Ce 4f states were treated in the open- 
core approximation and thus are characterized by integer number occupation 4 f 1 . We háve 
carefully checked that the calculation with these parameters converge. The hrst principles 
crystal ficld (CF) calculations were performed using the rnethod described in Refs 22 . Within 
this rnethod, the electronic structure and ground states charge density is obtained using the 
full potential APW + lo rnethod. The CF parameters of the microscopic CF hamiltonian 
originate frorn the aspherical part of the total single particlc DFT potential in the crystal. 
To eliminate the šelf interaction, the self-consistent proceduře is hrst converge with the 4/ 
electrons in core 22 , which is the open-core approximation ušed in this work. 

To obtain the dynamical properties the Hellmann-Feynman (HF) forces acting on the 
atoms starting from the equlibrium are required. For that we háve utilized the unspin- 
polarized density functional theory rnethod implementing the projector-augmented wave 
(PAW) formalism 23,24 to describe the electron-ion interactions. The approximation for the 
exchange-correlation potential was ušed the same as in WIEN2k calculations, i.e. GGA. 
The Kohn-Sham wave functions are expanded into plane waves up to a cutoff energy of 
600 eV. The valence electrons of cerium, lanthanum, palladium, and aluminium atoms are 
represented by conňgurations of 5s 2 5p 6 5d 1 6s 2 , 5s 2 5p 6 4/ 1 5d 1 6s 2 , 5s 1 4d 9 , and 3s 2 3p 1 , respec¬ 
tively. The calculations for the geometry optimazation ušed the 13 x 13 x 8(18 x 18 x 8) 
Monkhorst-Pack mesh of k-points for the orthorhombic(tetragonal) phase and the conver- 
gence eriteria for the systenťs total energy and residual HF forces were of 10 -7 eV and 10~ 4 
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eV/Á, respectively. 

The dynamical properties of the both LaPd 2 Al 2 and CePd 2 Al 2 lattice are obtained within 
the quasi-harmonic approximation and the direct method 25,26 , which utilizes the DFT cal¬ 
culated HF forces acting on all atoms in a given supercell. Phonons are calculated for the 
supercells composed of 80 atoms (16 Ce(La), 32 Pd and 32 AI atoms), which are created 
from the optimized unit cells. Such supercells were found to be large enough to avoid con- 
tributions from atoms bclonging to the periodic images, as confirmed by the clcments of the 
force constant matrices which decay by more than three orders of magnitude at the distances 
smaller than the boundaries of the supercells. The Brillouin zone integration is performed 
with the reduced number of /c-points 7 x 7 x 3(9 x 9 x 3) for the tetragonal (orthorhombic) 
phases. The non-vanishing HF forces required to construct respective dynamical matrix 
D(k) are generated by displacing the symmetry non-equivalent Ce(La), Pd and AI atoms 
from their equilibrium positions by the amplitudě of ±0.01 Á. Hence, for each configuration 
the total number of the calculated displacements amounts to 10 for both phases. 

III. RESULTS AND DISCUSSION 

The crystal field parameters and corresponding energy levels háve been calculated assum- 
ing the two crystal structures of CePd 2 Al 2 - high-temperature tetragonal (CaBe 2 Ge 2 -type) 
and low-temperature orthorhombic Cmme. The structural parameters háve been taken from 
experiment.' The results are summarized in Table I both for tetragonal and orthorhombic 
Ce positions in the corresponding crystal structures. Both CF hamiltonians are dominated 
by second order term B20. This is the reason that the full splitting of the Ce 3+ ground statě 
multiplet remains similar for both tetragonal and orthorhombic phases. In the tetragonal 
structure, the second excited statě lics at 11.8 meV what perfectly fits to energy of the hy- 
pothetical crystal field level assumed around 12 meV when analyzing the neutron scattering 
data 7 . The energy of this level changes only marginally when going from the tetragonal 
to the orhorhombic phase. The hrst excited level is calculated to be at 3.6 meV what is 
noticeably higher than the experimental value of 1.4 meV.' Possible explanation can be 
the complete neglect of hybridization in our hrst principles CF calculations. In compounds 
the hybridization of 4/ states with valence states is very sensitive to the details of position 
of the 4/ levels 2 '. This crystal field level then considerably shifts to higher energies when 
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B20 

B22 

B40 

B42 

B44 

Ai 

A 2 

tetragonal 

0.655 

0.0 

0.109xl0 -2 

0.0 

0.173x 10 _1 

3.6 

11.8 

orthorhombic 

0.628 

-0.389 

0.104xl0~ 2 

-0.405x 10“ : 

2 -0.127xl0 _1 

4.8 

12.2 


TABLE I: The calculated crystal field parameters for tetragonal and orthorhombic phase 
of CePd 2 Al 2 in SteveiTs notation Blin*0lm. The meV units are ušed for both Blm and the 

energy splitting A. 


CePd 2 Al 2 

a(Á) 

6(Á) 

c(Á) 

ZCe 

ZPd 

Z AI 

tetragonal 

4.435(4.415) 

- 

10.090(9.874) 

0.752 

0.372 

0.123 

orthorhombic 

6.274(6.268) 6.271(6.132) 10.089 (9.886) 0.753(0.756) 0.372 (0.371) 0.124(0.119) 

LaPd2Al2 

a(Á) 

6(Á) 

c(Á) 

ZLa 

ZPd 

Z AI 

tetragonal 

4.451(4.441) 

- 

10.087(9.896) 

0.753(0.749) 

0.371(0.376) 

0.121(0.138) 

orthorhombic 

6.334(6.313) 6.242(6.194) 10.105 (9.921) 0.752(0.745) 

0.371(0.377) 

0.129(0.149) 


TABLE II: The calculated lattice paramters a, b, c in Á and internal degrees of freedom 
(zx, X=Ce/La,Pd,Al). The experimental parameters in braekets for the Ce-phase and 
La-phase are taken from Ref. 7,19 and Ref. 19 , respectively. 


undergoing transition to the orthorhombic phase. This is in a qualitative agreement with 
the experimental observation'. 

Prior to the investigation of the lattice dynamics the geometrical optimizations need to 
be performed for the tetragonal PA/nmm and orthorhombic Cmme phases of CePd 2 Al 2 
and LaPd 2 Al 2 . Lattice parameters as well as the optimized internal degrees of freedom are 
summarized and compared to the experimental data in Table II. Calculated lattice constants 
agree well with the recorded ones, being slightly higher, but no more than 2%, than measured 
ones. This is a typical feature of the PBE functional for the exchange-correlations effects 
that tends to underestimate the binding energy. The overall agreement is therefore very 
good. Also the relaxed internal degrees of freedom, z-component of the Wyckoff positions 
(for details see Ref.') show excellcnt agreement. 

Lattice dynamics was investigated for both CePd 2 Al 2 and LaPd 2 Al 2 . The dispersion 
relations of orthorhombic and tetragonal structures are shown in Figs. 1,2, respectively. In 
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both structures there are no imaginary modes for any q-wave vector within the Brillouin 
zone and therefore the lattice structures are stable. CePd 2 Al 2 in both structures has higher 
slopes, see e.g. the ones between the Y and Z point, i.e along the [00((] direction that indicates 
higher clastic stiffness along the c-axis. Another difference concerns the degeneracy of the 
transversal modes in the čase of the tetragonal structure that is much higher in the čase of 
LaPd 2 Al 2 with respect to CePd 2 Al 2 , see transversal modes e. g. along Y to X and Y to M, 
respectively. 

The total and partial phonon density of states (DOS), for deňnition see e.g. Ref. 28 , 
for both CePd 2 Al 2 and LaPd 2 Al 2 are depicted in Figs. 3 and 4, respectively. The total 
phonon DOS for both compounds is very similar as regard frequency range owing to the fact 
that the rnass difference is also smáli (Mc e =140.115, Mdj= 138.906). Optical and acoustic 
frequencies are separated in both compounds and structures by a gap in frequency range 
between ca. 4.4-6.5 THz. The low frequency vibrations, in the range of 0-4.4 THz are 
dominated by the motions of the rare earth metal as well as palladium. The AI atoms 
vibrate at much higher frequencies. Here, the shape of the spectra for tetragonal structure 
is also very similar between both rare earths based compounds with higher frequencies of 
the CePd 2 Al 2 . In contrast the orthorhombic structure shows a large shift of phonon DOS 
at the higher frequencies (AI vibrations) for the Ce- vs. LaPd 2 Al 2 compound. Please notě 
the distinction of frequency range of vibrations for the AI atoms, i.e. a clear difference 
between motions of atoms with Wyckoff positions 2c and 2b and 4 b and 4 g of tetragonal and 
orhorhombic structure, respectively. Also, for the orthorhombic structure the center of rnass 
of the lower frequency of phonon DOS (R+Pd atoms) is shifted by ca 0.5 THz to higher 
frequencies for the CePd 2 Al 2 with respect to the LaPd 2 Al 2 . There is also a remarkable shift 
of the low frequency DOS (R+Pd atoms) to higher energies in the orthorhombic structure 
of CePd 2 Al 2 with respect to the tetragonal one. No such shift occurs for LaPd 2 Al 2 . 

The total phonon density of states could be ušed to estimate the Debye temperature (Od)- 
Using lower frequencies 0-2THz and fitting function of 0p(w) = acc 2 , where a is a parameter 
for the fitting function of©D = (yy) 1 / 3 , where N is number of atoms in unit cell. The total 
DOS of LaPd 2 Al 2 in both structures is similar, see Figs. 3 and 4, the estimated Debye 
temperatures are practically the same, Od = 185(183)/! for the orhorhombic (tetragonal) 
structure, respectively. Slighlty larger difference is in the čase of CePd 2 Al 2 , here O D = 189/í 
vs. 205 K between the tetragonal and orthorhombic one. Calculated Debye temperature for 
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LaPd 2 Al 2 is considerably larger than the one obtained from measurements of specific heat 
and elcctric resistivity, ( 0 £> — 117 — 127 K 29 , however in a very good agreement with the 
value obtained by Kitagawa et al., <d D = 174 K 18 . In fact, results of the models ušed in 
both references are in agreement considering that the former one considers only the three 
acoustic branches for the 0 o evaluation at very low temperatures while 0^ in the latter 
one corresponds to the total number of 15 phonon branches.Our calculation is in line with 
the latter approach. 

The calculated phonon DOS enable us to evaluate some thermodynamical quantities that 
are phonon-dependent. Here, we restrict our investigations to phonons calculated within 
the harmonie approximation (quasi-harmonic approximations) and analýze mean-squared 
amplitudes of atomic vibrations ( Uij ) ands the lattice contribution to the heat capacity as 
a function of temperature. Mean-squared vibrational amplitudes of a given atom in crystal 
lattice constitute a second-rank symmetric tensor and they can be expressed by the diagonál 
and off-diagonal partial phonon DOS 28,30 . In Fig. ?? the MSD in A 2 is clisplayed for the 
LaPd 2 Al 2 .. 

The calculated lattice heat capacity at constant pressure ( C p ) is compared with the mea- 
sured data in Fig. 6 for LaPd 2 Al 2 . The speciňc heat of LaPd 2 Al 2 is reproduced with a very 
good agreement including the low-temperature part (just above the T c ). The direct compar- 
ison for CePd 2 Al 2 is not possible because the specific heat of the Ce materiál includes also 
signiňcant magnetic contribution and a large heat capacity connected with the structural 
phase transition at 13 K. Neverthclless, the calculated heat capacity at constant pressure 
shows the negligible influence of the electronic contribution, to be linearly dependent along 
the temperature and corresponding to the number of electronic States at Fcrmi level (Som- 
merfeld coefhcient). This becorne more signiňcat over T^XXXK, whereas at all determined 
temperatures the lattice contributions are dominating as expected, see Fig. XXX. 

Let us ňnally relate all the calculated results to the inelastic neutron scattering data 7 what 
was a primary goal of this work. The simultaneous occurrence of the high density of phonon 
States and the CF level at the same energy as the key preriquisite for the magnetoelastic 
coupling and formation of the vibron states. Our present calculations are fully in agreement 
with such a scheme. The high phonon density of states around 12 meV is i) well in agreement 
with the experimental observation 7 and ii) appears at energy which is the same or very 
close to the calculated energy of the seconcl CF excited statě. We notě that there is no 
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CF excitation experimentally observed around 12 meV, but it is just the energy of the 
hypothetical CF level assumed when analyzing the neutron scattering data 7 . 

When undergoing the structural transition from tetragonal to orthorhombic structure 
at 13 K, the experiment revealed a strong shift of the hrst magnetic excitation to higher 
energies 7 . Such a shift is qualitatively well reproduced by our ňrst-principles electronic 
structure calculations although there is a certain quantitative clifference between the experi- 
mental and calcnlated energy as discussed above. The experimentally observed development 
of the energy of the hrst excited CF level can be thus primarily ascribed to CF changes when 
transforming from tetragonal to orthorhombic structure. 

Our calculations show also sorne clear changes of the phonon spectra when comparing 
the tetragonal and orthorhombic čase. As concern the low-energy part of the of CePd 2 Al 2 
phonon spectra, there is a clear shift of the DOS (R+Pd atorns) to higher energies in the 
orthorhombic structure with respect to the tetragonal structure. No such shift occurs for 
LaPd 2 Al 2 . This can point either to a different characteristics of the structural transition in 
CePd 2 Al 2 and LaPd 2 Al 2 (see also 19 ). 

The calcnlated lattice heat capacity at constant pressure (C p ) is compared with the mea- 
sured data in Fig. 6 for LaPd 2 Al 2 . The specific heat of LaPd 2 Al 2 is reprodnced with 
a reasonably good agreement. Rather snrprisingly, the calculation corresponding to the 
tetragonal lattice agrees almost perfectly with the experimental data whilc the data for the 
orthorhombic structure slightly exceed the measured values in the low-temperature part 
(just above the T c ). We remind that LaPd 2 Al 2 undergoes the structural transition to the 
orthorhombic phase around 90 K 6,12,19 , so we would expect better agreement with the or¬ 
thorhombic lattice calculation at low temperatures. The clirect comparison for CePd 2 Al 2 is 
not possiblc because the specific heat of the Ce materiál includes also signiňcant magnetic 
contribution and a large heat capacity connected with the structural phase transition at 13 
K. 

MSD, Cp, enegy difference with phonons and force constant, irreps 

IV. CONCLUSIONS 

The first-principles electronic structure calculations for CePd 2 Al 2 lead to crystal ficld 
excitations at 3.6 meV and 11.8 meV for the tetragonal structure and at 4.8 meV and 



12.2 meV for the orthorhombic structure. The change for the first excited levcl agrees 
qualitatively well with the experiment. The energy of the second excited lcvel fits perfectly 
with the model ušed to describe the inelastic neutron scattering data. 

The lattice dynamics calculations show that in both tetragonal and orthorhombic struc- 
tures there are no imaginary modes for any q-wave vector within the Brillouin zone and 
therefore the lattice structures are stable. There is a large phonon DOS at around 12 meV 
in agreement with the experimental observation. The calculation indicate a shift of the low 
frequency DOS (R+Pd atoms) to higher energies in the orthorhombic structure of CePd 2 Al 2 
with respect to the tetragonal structure. No such shift occurs for LaPd 2 Al 2 . 
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FIG. 1: Phonon dispersion relation of LaPd 2 Al 2 and CePd 2 AF in ortorhombic structure 

(Cinme, 67). 
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FIG. 2: Phonon dispersion relation of LaPd 2 Al 2 and CePd 2 Al 2 in tetragonal structure 

(PA/mmm, 129). 
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FIG. 3: Total and partial density of states in (states/THz) of LaPd 2 Al 2 and CePd 2 Al 2 in 

ortorhombic structure (Cmme, 67). 
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FIG. 4: Total and partial density of states in (states/THz) of LaPd 2 Al 2 and CePd 2 Al 2 in 

tetragonal structure (P4/mmm, 129). 
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FIG. 5: Mean square displacement of the LaPc^A^. The solid and dashed line mark for 
the orthorhombic and tetragonal structure, respetively. 
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FIG. 6: Comparison of measured and calculated speciíic heat of LaPd 2 Al 2 . The calculated 
curve includes, beside the dominant phonon contribution, also the electronic part taking 
the experimentally determined gamma coefficient of 15mJmol—1K—2 29 . 


17 







